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Abstract 

We compute the angular power spectrum of temperature anisotropics on 
the microwave sky in the cosmic texture theory, with standard recombination 
assumed. The spectrum shows 'Doppler' peaks analogous to those in scenarios 
based on primordial adiabatic fluctuations such as 'standard CDM', but at 
quite different angular scales. There appear to be excellent prospects for using 
this as a discriminant between inflationary and cosmic defect theories. 

The cosmic microwave background (CMB) anisotropy is the cleanest probe we have of 
structure formation in the universe, providing our best hope for definitively confirming or ex- 
cluding theories of the origin of structure. The COBE observation [1] provided the first solid 
evidence of anisotropy, but due to its limited angular resolution it provided little more than 
an overall normalisation and rough evidence of scale invariance. Higher resolution experi- 
ments offer the far richer prospect of a detailed measurement of the angular power spectrum 
and statistical properties of the anisotropics, which will give very powerful constraints on 
theories. 

Existing theories of the origin of structure fall into two categories, based upon either the 
amplification of quantum fluctuations in a scalar field during inflation or upon the process 
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of symmetry breaking and field ordering. The latter category includes the theories of cosmic 
strings, global monopoles and textures. While techniques for computing density perturba- 
tions in inflationary theories are well established [2,3], greater uncertainty is attached to the 
field ordering theories because they are inherently nonlinear and nonGaussian. 

Nevertheless there has been progress, particularly in those theories for which the non- 
linear sigma model describes the field ordering process [4]. Recently the degree scale CMB 
anisotropics were calculated for the cosmic defect theories under the simplifying assumption 
that the universe was fully ionised [5]. Reionisation is more likely in defect theories than in 
Gaussian theories because large density perturbations around the defects could cause early 
star formation, releasing ionising radiation. But the extent of this effect is uncertain because 
it depends on the efficiency of star formation, which is poorly understood. 

In this Letter we compute the power spectrum of CMB anisotropics in the texture theory 
assuming instead standard recombination. In particular, we wish to see whether acoustic 
oscillations of the photon-baryon-electron (PBE) fluid produce 'Doppler peaks' analogous 
to those in inflationary theories. The calculation is harder than in the reionised case, but 
simplifies if we concentrate on the power spectrum alone rather than attempting to make 
sky maps. While the nonGaussianity of such maps is undoubtedly a key signature for cosmic 
defects, the CMB anisotropy power spectrum is likely to be measured first and that is what 
we focus on here. 

The power spectra of perturbations depend only on the two-point correlations of the field 
stress energy tensor, 0^, and these one can hope to model in a simple way. We use the 
three dimensional field simulations to construct a model for the correlations, which are 
then fed into linear perturbation theory codes to compute density perturbations and CMB 
anisotropy spectra. As a check of our model, we perform full 3d fluid simulations of the 
source fields, gravity and matter in the tight coupling epoch and compare the perturbation 
spectra and cross correlations in the simulations with those predicted by the model. 

There are a number of advantages to this hybrid approach. A single 3d simulation has 
very limited dynamic range, but many such simulations can be combined to help construct 
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a model valid over all relevant length scales. Such a model allows one to calculate ensemble 
averaged power spectra without the noise inherent in 3d simulations. It can also be used in a 
Boltzmann calculation to compute the anisotropics for arbitrary ionisation histories. Finally, 
the model may provide some insight into the physics behind the generation of perturbations. 

In field ordering theories, the simplest assumption is that the universe began in a ho- 
mogeneous and isotropic initial state. A symmetry breaking phase transition occurs, during 
which some field with 'angular' degrees of freedom is given a nonzero vacuum expectation 
value, <f> 2 = <f>Q. The source for perturbations is 0^, the stress energy tensor of the ordering 
fields, interacting with the matter and radiation only gravitationally. Since the field ordering 
process is causal, it follows that the correlations of the fluctuating part of are strictly 
zero for spacetime points whose past light cones extrapolated back to the time of the phase 
transition do not overlap. 

We work in synchronous gauge, setting all perturbation variables zero before the sym- 
metry breaking phase transition. The metric and matter perturbation variables hij and 6^ 
obey causal evolution equations. Since the power spectrum is the Fourier transform of the 
correlation function, it follows that all these variables have 'white noise' power spectra on 
superhorizon scales at all times. 

The source fields interact gravitationally, as described by the linearised Einstein equa- 
tions. These are simplified by decomposing them into scalar, vector and tensor pieces: in 
Fourier symmetric tensor, Tij(x) = EkTf,-(k)e , can be written as 



where Tf hi = h{Tj 3 = T?kj = = 0. A complete set of evolution equations for the metric 
perturbation variables is: 
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hi + 2-hl + k 2 h% = 167TG0J. (4) 

hY + 2 a -hJ = 16.GeJ ■ (5) 
a 

where II = <9;©o;. The metric perturbations are defined by ds 2 = a 2 (r)(-G?r 2 + (Sij + 
hij(x,T))dx l dxi), with r conformal time and a(r) the scale factor, and k 2 h~ = k 2 (h — h s ) 
is proportional to the scalar curvature of the constant r spatial slices. The variables p^ } 
Pn } Sn and ujv refer as usual to the density, pressure, fractional density perturbation and 
velocity of each fluid - photons, baryons, cold dark matter and neutrinos. In the tightly 
coupled regime, the equations for the fractional perturbation in the photon density 8 1 are 

S\ = — k H — ikv 1 ; v 1 = -ikc 2 S 1 (1 — 3c 2 )u 7 (6) 

3 3 4 a 

where c 2 = |/0 7 /(/0 7 + Pb) is the square of the speed of sound, pb the baryon density. We 
model the neutrinos fluid with viscosity. 

In the 'stiff' approximation [6], is covariantly conserved with respect to the back- 
ground metric: 

/ 2 

©oo + -(0oo + 0) = n; II + 2-II = --(0 + 20 5 ). (7) 
a a 6 

We would like to model the behaviour of 0^ while respecting these constraints. The correct 

underlying dynamics is that of scalar fields obeying the nonlinear sigma model equation, 

but a crude description is that spatial gradients in the scalar fields are 'frozen' outside the 

horizon, and redshift away after horizon crossing. This may be viewed as a change in the 

effective equation of state: if there were only spatial gradients outside the horizon, we would 

have = —0oo- If these gradients redshifted away like radiation we would have O o = +0 

inside the horizon. This rough picture forms the basis of our model. 

In our model we treat the source as a fluid with a scale-dependent equation of state: 

= j(k, r)0 O o and 5 = 7 s ' (k, r)0 O o, where 7 and 7^ are determined from field simulations. 

In the pure matter or radiation epochs they have scaling form, being functions only of kr. 

With 7 and 7^ fixed, equations (7) determine the evolution of all scalar parts of 0^. 
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The main simplification this model makes is that the evolution equations for 0^ are 
linear. While a given mode is outside the horizon it quickly settles into a 'scaling' solution 
specified by a single amplitude. It follows that in the long time limit unequal time cross 
correlations factorise, i.e. 

(A k (r)£_ k (r')) = (A_ k (r)£ k (r')) = a(k, r)b(k, r') (8) 

where a(k, r) and b(k, r) are 'master' functions satisfying the same linear equation that the 
random field modes A k and B^ do. Setting B^ = A k and r = r', one sees that a 2 (&,r) is 
just the power spectrum of A k (r). The functions 7 and 7^ are given by (©oo0)/(©oo) anc ^ 
(0oo0 5 ) / (©oo) 1 which we measure in field simulations and perform a simple fit [7] to specify 
the model. 

The initial conditions for the model are determined by scaling and causality. Scaling 
implies that all correlators should be specified in terms of a single scale r alone. Dimensional 
analysis then gives 

(|0oo(k,r)| 2 ) = ^M (9) 

where V is a fiducial comoving volume. Causality implies a white noise spectrum for O o, 
so that /(fa) = constant for fa < 1. Thus for each k the master function for O o should be 
by a fixed constant times r~2 at early times. Scaling of and equation (7) fix 7(0) = — | 
in the radiation era. 

Let us make two comments on the most obvious limitations of the model. Examination 
of the large N expressions [8] shows that the Fourier modes of the energy momentum tensor 
far outside the horizon (fa <C 1), or far inside the horizon (fa ^> 1) are dominated by 
interference terms involving horizon wavelength modes, where the most important dynamics 
is taking place. These interference effects are unlikely to be well represented by the dynamics 
of a fluid in which all Fourier modes decouple. However, the dominant perturbations in a 
given Fourier mode of the metric and matter variables are produced by the source around 
horizon crossing, so all we really require is that the model adequately represents the source 
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at this time. At horizon crossing, only dynamics on a single length scale are involved, and 
these may be reasonably well described by fluid equations. 

A second concern regards unequal time correlations. Since the perturbation equations 
are linear, the solution for the set of perturbation vartiables 8 a takes the form 8 a (k } r) = 
J T dT'G a b(k } r, T')Sb(k, t'), where Sb are the sources O o, II , and G a b the appropriate Greens 
functions. It follows that all power spectra of the perturbations are determined by the 
unequal time source correlators. Equation (8) shows that the model actually builds in 
maximal unequal time correlations, leading to an overestimate of the coherence in time of 
the source terms [9]. We argue again, however, that around horizon crossing there is only one 
length and timescale involved, and the timescale for decoherence of the source is of the same 
order as that for the source to redshift away. This is what our model builds in. Nevertheless 
it is of some importance to check the model predictions against a full 3d simulation to see 
whether the predicted coherent oscillations are really there. 

Figure 1 shows the results of a 3d simulation of equations (2 -5) sourced by texture fields. 
The Figure shows the power in the variables h~ , <5 7 , 8c, u 7 from the 3d code against those 
predicted by model, at the instant of matter-radiation decoupling. There are no adjustable 
parameters in the comparison, and there is good qualitative agreement between the simu- 
lations and the model. We have similarly checked the power spectra of the source terms 
0oo and II and found reasonable agreement. As mentioned, we are particularly concerned 
to check whether the radiation oscillations evident in the model are really present, since 
these give rise to the Doppler peaks. The power spectrum for 8 1 shows some evidence of 
coherent oscillations, but of small amplitude. The power spectrum of u 7 shows larger os- 
cillations. The difference is explained by looking at the behaviour of the model. Whereas 
the master function for 8 1 oscillates about zero, with neighbouring peaks of similar ampli- 
tude, u 7 oscillates about a constant offset, so that neighbouring peaks have quite different 
amplitude. Smearing effects, whether from numerical effects or from real incoherence in the 
perturbations, do less to reduce the amplitude of these oscillations. 

A difficulty in computing power spectra from the 3d simulations is that oscillations in 
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FIG. 1. Power spectra from three dimensional simulations (solid lines) compared to those 
predicted by the equation of state model (dashed lines). The perturbation variables are those 
defined in equations (2-6). The vertical scales are arbitrary: large ticks are separated by an order 
of magnitude. 
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FIG. 2. The cross correlation k 3 (S 1 (k)h~(k)) measured in a single simulation (solid) is com- 
pared to the predictions of our model (dashed line). This provides strong evidence that coherent 
oscillations exist in the photon baryon fluid. 
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power are smeared out by sampling effects and by noise, which prevent one ever seeing a zero 
in a positive definite quantity like the power. A cleaner check of whether phase-coherent 
oscillations are present in the 3d simulations is to measure the cross correlation of 8 1 with 
some variable that does not have oscillatory behavior, such as h~ . Figure 2 shows (S*yh~) 
measured in the simulations compared with the model predictions. If the oscillations were 
incoherent this cross correlation would average to zero. Similarly, we have checked (6c h~) 
and confirmed that a single sign change takes place just as the model predicts. 

All this we take as evidence that the scalar part of the model is working well. The 
large scale microwave anisotropics also involve vector and tensor perturbations and we have 
treated these in the same 'factorisation' approximation. That is, we compute the power 
spectrum for Q v and T in 3d simulations and approximate the unequal time correlation 
functions as products of the square roots of the relevant power spectra. As mentioned 
above, it is hard to see oscillations in the true power spectra from those measured in the 
3d simulations. To check our sensitivity to this effect we have also modelled the vector and 
tensor cross correlations using oscillatory 'master functions' whose envelope has the same 
amplitude. The change alters the vector and tensor power spectra somewhat, but has a 
small effect on our final multipole power spectrum. 

In the approximation of instantaneous recombination one can write the full expression for 
the microwave anisotropy in a direction n as the sum of 'intrinsic', 'Doppler' and 'integrated 
Sachs- Wolfe' pieces, 

y-(n) = ^a; m YJ m (0, <f>) = - v 7 • n - ^ J drh^n'ri 1 (10) 

Converting this into a spectrum of C[s involves ensemble averaging which is simple in the 
factorisation approximation. The scalar expression is well known (see e.g. [10]), the vec- 
tor and tensor contributions are even simpler. We define the vector master function by 
(hy (k)hj(-k)) = h v '(k,T)h v '(k,T')(6 tJ -k t kj), and the tensor function by (A^(k)A^(-k)) = 
h T (k, r)h T (k, T r )(8ikS J i + ...) where ... denotes terms imposed by the transverse-traceless and 
symmetry conditions. Some algebra then shows that 
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<i«Li 2 > = 16*7 ^ ,{i+i}{ rt 3 ' {kT)hv{k - T}f (u) 

(l-Ll^^/^illl,/^,,,,,^,^), 2 (12) 

The latter equation agrees with (and simplifies) the formula given by Abbot and Wise 
[11]. The functions h v and h T are obtained by solving the Einstein equations (4) and (5) 
numerically, using the master functions for Q v and T inferred from 3d simulations. 

The final CMB anisotropy power spectrum is presented in Figure 3, along with various 
components. The results shown were computed using a fluid code, in the 'instantaneous 
recombination' approximation. We have checked the scalar and tensor calculations using a 
full Boltzmann code, with very good agreement. The cosmological parameters used were 

h = 0.5, o = i, n B = 0.05. 

A word about uncertainties. We have found the shape of the spectrum for / > 100 to be 
remarkably insensitive to the details of our model: plausible variations change the results 
by less than 10%. Likewise the qualitative shape of the vector and tensor contibutions is 
fairly model-independent. We attach greater uncertainty (of order 50%) to amplitude of 
the perturbations on large (COBE) scales, first because of the complexity of extracting the 
vector and tensor master functions from the 3d code, discussed above, and second because 
the scalar power on these scales is sensitive to variations in the 'equations of state' as one 
passes from the radiation to the matter era. More detailed numerical studies will resolve 
these uncertainties. 

The most striking difference between the texture and inflationary spectra is that 
'Doppler' peaks and troughs are almost 'out of phase' - where inflation predicts a maxi- 
mum texture gives a minimum, and vice versa. This behaviour is reminiscent of that found 
in 'isocurvature' models [12], and it occurs for a similar reason [13]. The radiation oscilla- 
tions are driven by the metric perturbations: from equations (6) we have S 1 + c|<5 7 = — \h. 
The phase of the oscillations is determined by the behaviour of h as a mode crosses the 
horizon. In 'adiabatic' theories like the simplest inflationary models, one has superhorizon 
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FIG. 3. The predicted anisotropy power spectrum C\ = (|a/ m | 2 ) (equation (11)), decomposed 
into scalar, vector and tensor contributions. The vector and tensor parts have been multiplied by a 
factor of two for clarity. Shown for comparison (grey curve) is the prediction of the simplest infla- 
tionary theory with cold dark matter with the same cosmological parameters, h = 0.5, 12/; = 0.05, 
and = 1. 
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curvature perturbations, in which the metric perturbation h oc r 2 . Thus the force driving for 
the radiation oscillations is to a first approximation constant. In the 'isocurvature' theories, 
the source term is not constant and this leads to a relative phase shift in the oscillations. 

The field ordering theories are both adiabatic and isocurvature: there are no physical 
perturbations in the space curvature or entropy on superhorizon scales. This leads to a 
different behavior of the h forcing term around horizon crossing. Scale invariance and 
dimensional analysis imply that 



with g a dimensionless function. In standard 'adiabatic' theories, h oc r 2 so g oc kr for 
small kr. But for field ordering theories, the small k behaviour of g is fixed instead by 
causality: h must have a 'white noise' power spectrum, and so g = constant and thus 

3 1 

h oc T2 at small kr. Thus the forcing term decreases as r~2. It follows that S 1 reaches 
its first maximum sooner than in the adiabatic theories. It is intriguing that, at least in 
the context of adiabatic theories, the large scale angular power spectrum, and in particular 
the phase of the 'Doppler' peaks could be telling us something as fundamental as whether 
the perturbations were generated causally within the standard big bang, or were necessarily 
generated in a preceding inflationary epoch. 

We thank D. Coulson for collaboration in development of the 3d code. This work was 
partially supported by NSF contract PHY90-21984, and the David and Lucile Packard Foun- 
dation. 
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